Remarks about the thermodynamic hmit in selfgravitating systems 
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The present effort addresses the question about the existence of a well-defined thermodynamic 
limit for the astrophysical systems with the following power law form: to tend the number of 
particles, A'', the total energy, E, and the characteristic linear dimension of the system, L , to 
infinity, keeping constant E /N^^ and L/N^^ , being A_e and A_t certain scaling exponent constant. 
This study is carried out for a system constituted by a non-rotating fluid under the influence of 
its own Newtonian gravitational interaction. The analysis yields that a thermodynamic limit of 
Qf) • the above form will only appear when the local pressure depends on the energy density of fluid as 

' p = 7£, being 7 certain constant. Therefore, a thermodynamic limit with a power law form can 

, be only satisfled by a reduced set of models, such as the selfgravitating gas of fermions and the 

^Nj ■ Antonov isothermal model. 
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\f-^ i In a series of recent papers, de Vega, Sanchez and Laliena discuss about the existence of a well-defined thermody- 
namic limit for the self-gravitating systems 0, 0, IE 13 • The first two authors claimed that the thermodynamic limit 
of a self-gravitating system can be taken by letting the number of particles, iV, and the volume, V, tend to infinity 
keeping the ratio N/V^ constant. On the other hand, Laliena exposes a series of reasons which in his viewpoint makes 
invalid the consideration of the above thermodynamic hmit. The present paper is aimed to contribute to clarify the 
question about the existence of a well-defined thermodynamic limit for self-gravitating systems. 

Our analysis starts from the consideration of a mean field approximation for the estimation of the microcanonical 
volume T4^ of a non-rotating fluid under the infiuence of its own Newtonian gravitational interaction: 

o ; 
cn . 

O ; Wmf {E, N,L)^ / Vp (r) Ve (r) V<j> (r) S (E - H[e, p, 0]) (5 (iV - N[p]) 5 {<f> (r) - G [p, r]} exp {^[e, p]} , (1) 

^ ■ 

^p^, being p (r) and e (r), the particles and the internal energy density of tlie fluid at tire point r, and <j> (r), the Newtonian 
Q ■ potential at the same point. The functionals: 

: S[e,p]^ J d'r s{e{r),p{r)}, H[e,p,cj,]= J d'r e {r) + ^mp (r) cj, (r) , N [p] ^ J d^r p (r) , (2) 

• I— I . 

, represent the total entropy, energy and particles number for a given profile with e (r) and p (r) ; and the functional 
is the Green solution of the Poisson problem: 



A(j) = iTiGmp. (4) 

This approximation was analyzed in details in ref.j^. This system has been enclosed in a spherical rigid container 
with characteristic linear dimension L in order to avoid the long-range divergence of the Newtonian gravitational 
potential. On the other hand, the short-range divergence of this interaction has been also avoided by the consideration 
of quantum effects or the nature size of the particles which are compose the fluid, which are taken in an implicit manner 
in the entropy density of the fluid, s {e (r) , p (r)}. 
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Wmf {E, N, L) can be rewritten by using the Fourier representation of the delta functions: S {x) = f^^^ ^ exp {zx), 
being z — e + ik, yielding: 

Wmf {E, N)^ / — ^ / Vp (r) Ve (r) V4> (r) Vh (r) exp {C [e, p, 2, zi, J]} , (5) 

J-00 J-00 (^71") 

where z — P + ik and zi = /i + with 77 G M, being the functional £ [e, p, (j); z, zi, J] defined by: 

C[e, p,4>; z, z,, J] ^ S [e, p] + z {E ~ H [e, p,cf>]) + z,{N - N [p]) + J * i<f> ~ g [p]) . (6) 
The functional term 



J*{cf>-g[p]) = I J{r){c^{r)-g[p;r]} 



appears as consequence of the Fourier representation of the delta functional S {(p (r) — g [p; r]}. Here, J (r) is a complex 
function, J (r) = j (r) + ih (r), with j (r) G K.. 

The direct integration of the functional integrals in the expression ((SJ is a formidable task which could be only 
carried out by using adequate approximations. The usual methodology which is applied in this situations is the called 
steepest decent method. The application of this method is based on the asymptotic behavior of the thermodynamical 
variables and potentials in the many particle limit N ^ 1, which leads to estimate that the main contribution for 
the expectation values for the physical observables in the microcanonical ensemble comes from the most probably 
configurations of the system. Thus, the Boltzmann entropy of the system, Sb = In W , is obtained by using a min-max 
procedure: 



Sb [E, N, L) ~ min <^ max C [e, p, 0; /3, p, j] \ , (7) 

whose the stationary conditions leads to the structure equations of the equilibrium configurations. According to the 
definition of the functional C[e, p, (j); 13, p, j] , the equation © , the reader may recognize the Lagrange formalism of the 
Classical Thermodynamics. 

The presence of an additive kinetic part in the Hamiltonian of certain system leads to an exponential growing of 
the microcanonical volume W with the N increasing, and therefore, the Boltzmann entropy will grow proportional 
to N in the many particles limit, Sb = InVF oc TV. Therefore, let us now concentrate our attention in the analysis of 
the A'^-behavior of the thermodynamical variables with the growing of the number of particles. 

The usual thermodynamic limit for the extensive systems: 



E N 

N — + 00, keeping constant — and — , (8) 

where V is the volume of the system, is directly related with the extensive properties of these systems when the 
thermodynamical variables of the system are scaled by some scaling parameter a as follows: 



N ^N{a)= aN^ \ 

E^E{a)=^aE', \ ^ W W {a) = exp {a\nW) . (9) 
V ^V {a)^ aV, J 

In analogy with the extensive properties of the traditional systems, we will analyze the necessary conditions for the 
existence of the following power law self- similarity scaling behavior of the microcanonical variables E, N and L for 
the selfgravitating systems: 



E ^ E{a) ^ a^'^E, \ ^ W ^ W (a) ^ exp (alnW) , (10) 
L L {a) = a^^L, J 
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where A^; and A^, are certain constant scaling exponent which lead to an extensive character of the Boltzmann entropy. 
This kind of self-similarity behavior is directly related with a thermodynamic limit with a power law form: 



N ~f oo, keeping constant and . (11) 

The existence of this kind of self-similarity condition allows a considerable simplification of the thermodynamical 
description: the study can be performed by setting = 1 and considering the A^-dependence in the scaling laws by 
taking a = N. This scaling behavior is very useful in numerical experiments, since it allows us to extend the results 
of this kind of study on a finite system to much bigger systems. Contrary, the nontrivial A^-dependent behavior of the 
thermodynamical variables and potentials leads to a complication of any kind of study. 

In order to satisfy this scaling behavior for the global variables E, N, L and the Boltzmann entropy Sb, the local 
functions e (r), p (r), (p (r) and s (e, p; (j>) should be scaled as follows: 



) (j) (a) — a-^'>'( 

p{a) ^ a^'-p} ^ \ e^e{a)^a^^e }. (12) 
s — + s (a) — a^'^ s 



Since the characteristic particles density behaves as pc ^ N/L^, the scaling exponent for the particles density is 
Ap = 1 — 3Al. From the expression of the Newtonian potential (j2Jl is derived that its characteristic unit is (f)c ~ PcL^, 
and therefore, A^ = 1 — A^. The energy scaling exponent is equal to the scaling exponent of the total gravitational 
potential energy, so that, Ae = 2 — A^. The other scaling exponents are obtained by using identical reasonings. All 
these scaling exponents depend on the scaling exponent A^ as follows: 



Ap = 1 - 3Al = As, A0 = 1 - Al, 
A, = 2-4Ar, AB = 2-Ar. 



(13) 



In order to satisfy these scaling laws is also necessary that the entropy density exhibits to the following scaling 
behavior: 

s{a^'e,a^>'p) ^ a^>'s{e,p) , (14) 
This scaling property is satisfy if s (e, p) obeys to the following functional form: 



s(e,p) = pF(e/p"), (15) 

where ij = Ag/Ap, being F (x) an arbitrary function. This functional form leads to a simple relation between the local 
pressure p and the internal energy density e of the fluid. By taking into account the functional form (|15|) and using 
the relations 



d d d 

P= ttS (e, p), (3p = s (e, p) - e—s (e, p) - p—s (e, p) , (16) 
at ae Up 



it is straightforward derived the relation: 



p = 7e, (17) 

where 7 = 77 — 1. There are some well-known Hamiltonian systems which satisfy this kind of relation, as example, the 
system of nonrelativistic or ultrarelativistic noninteracting particles, without matter if they obey to the Boltzmann, 
Fermi-Dirac or Bose-Einstein Statistics. The scaling parameter A^ and A^ are obtained from the parameter 7 as 
follows: 



(18) 
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This result evidences that the power laws form for the self-similarity conditions (|10|l can be only satisfied by a 
reduced group of models whose microscopic picture obeys to the relation IjlTI) . Since 7 = § for the ideal gas of 
particles, the scaling exponents for the Antonov problem and the selfgravitating fermions model 0, ^| are given 
by — I and A^ = — -i, and therefore, they obey to the following thermodynamic limit: 



Cj 1 

iV — > cx), keeping constant — - and LN ^ . (19) 

This thermodynamic limit was established in the ref. @ for the self-gravitating nonrelativistic fermions by using other 
reasonings. Since the selfgravitating relativistic gas and classic hard sphere model [lllIT^ do not obey to the relation 
(|17|l . they do not posses a scaling behavior with a power law form (|10|l . 

This result clarifies a question about the discussion of the adequate thermodynamic limit for selfgravitating systems. 
All above selfgravitating models converge in the low density limit in the selfgravitating gas, and therefore, they exhibit 
a ^-growing of the energy with the 7/3 power law. However, they diverge in regard to the high density limit since 
them used different regularization procedure for the short-range divergence of the Newtonian potential. Generally 
speaking, such power law thermodynamic limit Hll|) does not exist in selfgravitating systems for the whole values of the 
microcanonical variables. Even the selfgravitating gas of fermions exhibits this behavior only in the nonrelativistic 
limit since the relation (|15|l disappears when the relativistic effects are taken into account for a massive enough 
selfgravitating system. 
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